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Abstract 

In this paper we propose a procedure which allows the construction of a large 
family of FIR d x d matrix wavelet filters by exploiting the one-to-one cor- 
respondence between QMF systems and orthogonal operators which commute 
with the shifts by two. A characterization of the class of filters of full rank type 
that can be obtained with such procedure is given. In particular, we restrict 
our attention to a special construction based on the representation of SO(2d) 
in terms of the elements of its Lie algebra. Explicit expressions for the filters in 
the case d = 2 are given, as a result of a local analysis of the parameterization 
obtained from perturbing the Haar system. 

Keywords: Full rank matrix filters, Multichannel wavelets, Quadrature Mirror 
Filters, Vector subdivision schemes 
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1. Introduction 

Parameterization of orthogonal and biorthogonal filters has been an impor- 
tant topic of research in the context of wavelet analysis (for the relation between 
filters and wavelets see e.g. |13j). Pioneering works in this area have been carried 
out, for example, by Pollen [TU] and Holschneider [7J, who proposed parameter- 
izations based on loop group factorization. Later Sweldens [12| introduced the 
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lifting scheme which is now widely used to construct biorthogonal families of fil- 
ters. For yet another interesting approach to the explicit construction of wavelet 
filters based on their correspondence with representations of the so-called Cuntz 
relations we refer the reader, for instance, to [8]. 

The aim of this paper is to extend the approach given in [7] to the construc- 
tion of an infinite family of orthogonal perfect reconstruction matrix filter banks 
of full rank type. These kinds of filters, introduced in [IHS], are associated to 
the multichannel wavelet analysis of signals which are vector- valued, consisting 
of several components typically associated to different channels. Signals of this 
type arise naturally in many application contexts; typical examples are: brain 
activity (EEG/MEG) data, colour images, multisensor data, financial time se- 
ries, etc. A classical scalar wavelet analysis applied the single components of 
this kind of data might not be appropriate, because it ignores the possible corre- 
lation among the channels. Multichannel wavelets provide a more effective tool, 
able to process such type of signals as "complete" vectors, so to extract peculiar 
information from the overall behaviour of the components and to reveal/exploit 
inter-correlations. 

The main challenge in the context of multichannel wavelet analysis is the 
construction of matrix filters satisfying the quadrature mirror filter (QMF) con- 
dition. In particular, the direct construction of the filter associated to the matrix 
scaling function from the QMF constraints is not convenient because of their 
non-linear nature. On the other hand, a Daubechies-like approach only gives 
rise to filters which are essentially diagonal, in the sense that they can be re- 
duced, via similarity transformations, to scalar schemes. In 4 a construction 
of orthogonal matrix refinable functions has been proposed as a spectral factor- 
ization problem, based on the connection between orthogonal and interpolatory 
matrix subdivision schemes. Nevertheless this approach cannot in principle give 
rise to any closed parameterized form of the filters. As to the construction of 
the corresponding wavelet filter, unlike the scalar case, there is no trivial al- 
ternating flip trick to express it in terms of the scaling function filter, due to 
the non-commutative nature of the matrix filter context. The numerical scheme 
presented in [5] is an example of effective numerical approach to the problem, 
which, nevertheless, relies on the a priori knowledge of the scaling filter. 

In this work, we address both the problems (construction of the scaling 
function and of the wavelet) at the same time, by deriving a procedure which 
provides a natural and intrinsic characterization of (quasi) all matrix QMF 
filters. The main idea is to exploit the one-to-one correspondence between 
matrix QMF filters and orthogonal operators that commute with translations 
by 2, so that starting with a trivial system and repeatedly applying such type 
of operators it is possible to produce large classes of filters. 

It is important to point out that the QMF filter construction in the ma- 
trix case can be viewed also in a multiwavelet perspective (see, for example, 
[51 US H])) but all the constructions proposed for this kind of bases do not 
take into account the full rank requirement (since multiwavelets still generate 
mult iresolut ion analyses for spaces of scalar rather than vector-valued func- 
tions). Our construction, on the other hand, can be straightforwardly applied 
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to the realization of multiwavelet bases. 

The paper is organized as follows. Section 2 presents some preliminaries and 
notation. In Section 3 we give a detailed exposition of our construction and a 
characterization result. Finally, in Section 4 we examine the particular case of 
dimension 2 providing explicit descriptions of families of parameterized filters, 
in particular those associated to the representation of the Lie group of rotations 
in terms of infinitesimal generators. 



2. Notation and basic facts 

By L 2 (Z) dxd , d e N, we denote the space of square-summable d x ci-matrix 
valued sequences, that is 

L 2 (Z) dxd = L:Z-> R dxd : \Mk)\l) < oo 

{ \fcez / 

where | • I2 denotes the 2-norm on the matrix space M <lxd . For notational sim- 
plicity we write L 2 (Z) d for L 2 (Z) dxl and L 2 (Z) for L 2 (Z) 1 . 

For two matrix sequences A,B e L 2 (Z) dxd , we define the inner product as 

(A,B)=J2Mk) T B(k) 

and we say that the two sequences are orthogonal to each other if 

(A,B)=0, 

where denotes the d x d matrix all of whose entries are zero. 

In particular, given two vector sequences a = ([ai(fc), . . . 7 a,d(k)] T : k € Z), 
b = ([&! (jfe), . . . , b d {k)] T : keZ)inL 2 (Z) d , then their inner (scalar) product is: 



< a, b >= a i( k )M k )- 



fcgZ i=l 

For n G Z, we define the Dirac delta sequence S n as the sequence whose 
elements satisfy: 

1, k = n 
0, k + n 

If n e Z, we introduce the translation and downsampling operators T„ and 
D acting on a sequence A € L 2 (Z) dxd respectively as: 



5 n (k) 



T n A = A(- — n), DA = A(2-). 

The transpose of the downsampling operator is the upsampling operator, de- 
scribed by the action 



(D T A)(k) 



A(k/2), k even 
0, k odd 
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The d x d identity matrix is denoted with J, and its columns, representing 
the canonical basis in M. d , with ej.,. . . , e<j. 

Two filters A and B are said to define a matrix QMF ( Quadrature Mirror 
Filter) system, if they satisfy the orthonormality conditions: 

A(jfA(j -2k)=J2 B(j) T B(j - 2k) = 25 (k)I, keZ, 

Y / A(j) T B(j-2k) = 0, keZ, 

which, using the inner product notation, may be written as: 

(A, T 2k A) = (B, T 2k B) = 25 (k)I, (A, T 2k B) = 0, k e Z. 

For a signal s € L?(Z d ) the application of a matrix valued filter is defined in 
the natural way through the convolution operation: 

A* s = M- - k) T s{k). 

feez 

A QMF system provides a perfect analysis/reconstruction scheme through the 
following two equations: 

s h- > (sq, Si) = (D(A * s), D(B * s)), D T (A* s ) + D T (B * s x ) = s. (1) 

Here we have used the notation A — A T {— ■). 

Instead of working with matrix valued functions we may also take the follow- 
ing equivalent picture by considering vector valued sequences. In analogy to the 
scalar case, a QMF filter system is generated by the dilates and translates of 2d 
vector- valued function a±, . . . , ad, and b\, . . . , bd, satisfying the orthonormality 
property: 

(a i ,a j (--k)) = 25 (k)S (i - j) 

(bi,bj(--k)) = 25 (k)5 (i- j), i,j = l,...,d,keZ 
(o<, &,•(•-*)) = 0. 

The collection of all these functions is then an orthonormal basis of L 2 (Z) d . We 
identify these vectors with the previously introduced matrix valued signal by 
taking them as their column vectors (analogue for b f-> B): 

d 

A = [a ly ...,a d ] = } j a j ef . 

3 = 1 

The decomposition and synthesis of a signal in L 2 (R d ) with respect to this 
orthonormal basis can then be obtained with the filters A and B as in (Til). 
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We conclude this section by introducing a crucial feature which is required to 
a matrix QMF system to be connected to a proper wavelet analysis context. A 
QMF system is said to satisfy the full rank condition if the sequence A satisfies: 

£A(2j) = £A(2j + l)=I. (2) 

This, in particular, implies that the coefficients of the sequence B sum up to 
the zero matrix. 

It can be shown, that full rank is a necessary (but not sufficient) condition 
for these filters to actually define a multichannel multiresolution analysis pQ, 
which is a natural extension of the well-known multiresolution analysis to vector- 
valued functions. Any (orthogonal) multichannel MRA is, in fact, generated 
by two square-integrable d x d-matrix valued functions F and G, representing, 
respectively, the matrix scaling function and the multichannel wavelet, satisfying 
the matrix two-scale relations: 

f = "£f(2- -j)a(j), g = J2f&- -j)b(j). 

In such context, the two-scale coefficient sequences A, B define a matrix QMF 
system and must possess the full rank property. In particular, this is related to 
the fact that the matrix wavelet G has vanishing zeroth moment: 

G(x) dx = 0. 

We remark that, as in the scalar case, further p—1 vanishing moments (usually 
required in applications) can be possessed by G, that is, 

/ x n G(x) dx — 0, n = l,...,p—l, 

only if the coefficients in the sequence A satisfy the additional sum rule condi- 
tions: 

]T(-l) fe fc"A(fc) = 0, n=l,...,p-l. 

feez 

3. The construction of matrix QMF systems 

The basic idea of the construction is easiest explained in the well known 
QMF setting of a scalar valued signal. 

Let a, b £ L 2 (Z) be such that the set of translates by 2Z: 

{o(- - 2m),b(- - 2m) : m £ Z} (3) 

is an orthonormal basis of L 2 (Z) , so that the two sequences form a QMF system. 
Note that the distinction between a and b is somehow arbitrary. However in 
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actual filtering applications one would like a to be the scaling function filter and 
b to be the wavelet filter. This would impose additional constraints (vanishing 
moments, etc) on these sequences. 

A trivial example of a QMF system in this purely algebraic setting would be 

a = Sq, b = Si. 

Consider now an orthogonal linear operator U acting in L 2 (Z) with the property 
that it commutes with the shifts by 2, that is: 

[U,T 2 ] := UT 2 - T 2 U = 0. 

Since U is orthogonal, the image of the basis ([3| is again a basis. Since U 
commutes with the shifts by 2, it has the following structure 

{{Ua)(- - 2m), ([/&)(• - 2m) : m € Z}. 

Therefore this basis is again a QMF system. 

Vice versa, any QMF system can be obtained from the trivial system in this 
way. Indeed, there is exactly one linear operator such that: 

fen i->- a(- - 2n), £ 2n+ i ^ &(' - 2n), 

since any linear operator is defined by its image of a basis. Therefore there 
is an explicit one-to-one correspondence between QMF filters and orthogonal 
operators that commute with translations by 2. Note that any such operator is 
defined through its image of So and <5i . The finite impulse response (FIR) filters 
correspond exactly to those operators which leave invariant the non closed linear 
subspace of compactly supported sequences. 

Therefore, if we can construct a manifold of such operators, we automatically 
have a family of QMF by applying them to any known QMF, in particular the 
trivial one {So, Si}. 

A family of such operators can be obtained in the following way (see Fig[T]). 
Pick an arbitrary 2x2 orthogonal matrix M <E 0(2). Let the image of Sq be 
its first column put at positions 0, 1 extended with zeros on both sides. Let the 
image of S\ be the second column, put at positions 0,1 and extended with zeros 
on both sides: 

So Mi^So + M 2 ,iSi, Si H> M\,%So + M 2 ^ 2 S X , 

Then extend mapping to a linear operator to all of L 2 (Z) by requiring that it 
commutes with the shift by two. This defines the operator U~o(M). 

To put it differently: for a sequence s = (s(n) : n e Z) e L 2 (Z), apply 
the same orthogonal matrix M to the vectors [s(2k), s(2k + 1)] T , k G Z. It 
is plain that this linear operator conserves the norm and that it satisfies the 
commutation property. 

Now, instead of taking this grouping, we may also consider the grouping 
[s(2k — l),s(2fc)] T . This defines a new operator, Ui(M), which can also be 
written as 

U x {M)=T_ x Uo{M)T x . (4) 
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Figure 1: Schematic view of the construction. The void circles are 0. The black dots corre- 
spond to d dimensional vectors. Note how the support grows with each application. 



The product of two such operators is again orthogonal and satisfies the com- 
mutation property. Thus, for any finite sequence of matrices M^ G 0(2), 
k = 0, 1, . . . K — 1, we have an operator 



[/(Mo, Mi, ... , M K -i) = U eK _, {Mk-!) ■ . . . ■ U ei (Mj) U €o (M ), (5) 



with eo, . . . , ck-i € {0, 1}, that can applied in order to obtain a parameteriza- 
tion of a family of QMF. 

In this case this parameterization is actually complete. Indeed, for a QMF 
of length < 2N the outermost blocks are all orthogonal to each other as follows 
from the orthonormality relations. Therefore, the spaces spanned by the left 
block of the a sequence and the left block of the b sequence are collinear. The 
same holds for the right blocks. Since they are orthogonal, there is a rotation 
such that the outermost coefficients can be sent to zero and thereby the length 
of the QMF sequence is reduced by 2. This gives an induction argument for 
proving the completeness of the parameterization. 

A further generalization is possible. Consider any fundamental domain F 
for the quotient group Z/2Z. By applying the same matrix operation to the 
grouping defined by F, we obtain again a linear operator that conserves the 
norm and commutes with the shifts by 2. Thus, for any sequence of fundamental 
domains and orthogonal matrices, we may consider the product 



It is clear that the same construction can be used to obtain all complex valued 
filters. It is enough to replace the orthogonal group and the orthogonal matrices 
by the unitary analogous. 

We now extend the construction of such operators to the vector case. As 
in the previous section, we may identify orthogonal operators in L 2 (Z) d that 
commute with T 2 and QMF systems. The trivial vector QMF system is given 



U FK {M K )-...-U Fl (M 1 )U Fa (M ). 




et6i, l = l,...,d}. 
We identify 4 submatrices 

" M M Mi, 2 " 
_ M 2 ,i M 2i2 



M 
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The orthogonal operator defined in L 2 (Z) d is acting on the sequences as follows: 

eiS Mi^ej^o + M 2 ,ie^i, ej5i ^ Mi^etfo + M 2;2 e^i. 

Again we extend this to all of L 2 (Z) d by requiring that the operator commutes 
with T 2 . 

In the language of blocks, this operator acts on a vector-valued signal s e 
L 2 (Z) d as follows. For each n e Z, we identify the two vectors s(2n) and 
s(2n + 1) with the 2<i-dimensional vector [s T (2n) s T (2n + 1)} T . The operation 
on the two vectors in the block is now multiplication with the matrix M. A 
second family is obtained by considering, for example, the grouping (2n — 1, 2n) 
which corresponds to another fundamental domain of Z/2Z. As before these 
operations may be combined to obtain a family of such operators and, a fortiori, 
a family of matrix filters. 

More precisely, in matrix notation, starting from a matrix QMF system 



{A(- - 2m), B(- - 2m) : m <= 
we apply a matrix M e 0(2d) to the sequence 



(6) 



or to the sequence 



A(2n) B{2n) 
A(2n+1) B(2n + l) 



A{2n-\) B(2n-l) 
A(2n) B(2n) 



and recursively do the same operation to the transformed sequences. 

It is convenient to introduce the description of the QMF systems in terms 
of symbols. For this, we associate with the matrix filter A the Laurent series 
matrix 

A(z) = Y J Mk)z\ zeC\{0}, 

which in fact is a Laurent polynomial matrix, if we assume that the filter has 
finite support. 

The so called subsymbols are defined through 



A l (z) = J2A(2k + l)z k , 1 = 0,1, 



and the following relation holds between the symbol and the subsymbols: 

A(z) = A (z 2 ) + zA 1 (z 2 ). 

Consider now also the symbol B(z) of the wavelet filter and the respective 
subsymbols Bq(z), B\{z) and define the following 2d x 2d polynomial matrix: 



L(z) = 



A (z) B {z) 
A^z) B^z) 
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which is also known as the polyphase matrix. 
Then the QMF conditions can be written as 

[L(z 2 )] H L(z 2 ) = 21, 

while the full rank requirement is equivalent to 



(7) 



L(l) 



I I 

I -I 



Straightforward computations show that the operator U(M) and the translation 
operator T\ are respectively acting on the symbols as follows: 



U(M) : L(z) i y ML(z), 



Ti : L(z) i-» 



z 



I 

- 1 ! 



L(z), 



(8) 
(9) 



so that the operator described in Q takes the following form in terms of sym- 
bols: 



zl 

1 



M 



I 

z- 1 ! 



" 


zl ' 


M K 


I " 


... M 2 


' 


zl ' 


Mi 


/ 


I 





z- 1 ! 







z~ x I 



Given a sequence of matrices Mk € 0(2d), k = 0,1,..., K, by alternately 
applying ^ and ([9| to the matrix L(z), we obtain the following family of QMF 
symbols: 



M L(z). 



We conclude the Section by giving an almost characterization of the set of 
matrix QMF that can be obtained using the above construction. 

Let J be the smallest interval that contains the support of all sequences in 
the QMF system. The even blocks correspond to grouping the elements with 
indices in {2n, 2n + 1} whereas the odd blocks are the one with {2n — l,2n}. 
The outermost blocks of an even or odd covering of Z are the left most and right 
most that intersect the support. The dimension of a block is defined to be the 
dimension of the vector space generated by the 2d vectors with indices in the 
block. 

Theorem 1. All QMF systems that can be obtained through this construction 
have the dimension of their outermost blocks < d. Vice versa all matrix QMF 
system of finite support length for which the dimension of outermost block is 
equal to d can be obtained through this construction. 

Proof. For length equal to 2 nothing has to be proved. 

That the dimensions of the outermost blocks are at most d is clear, since it 
is obtained by applying a linear map to a vector space of dimension at most d. 

Vice versa, suppose that the outermost blocks for either an even or odd 
covering are each d dimensional. By orthogonality of the basic vector sequences, 
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all the extremal blocks on the left are orthogonal to all the extremal blocks on 
the right. Since the dimension of each of these blocks is d it is possible to find 
an orthogonal matrix such that the left block is mapped into a block having 
only zeros to the left, and the right block is mapped into a blocking having 
zeros to the right. Therefore the length of the filters is reduced by 2. Now, the 
dimension of each of the new extremal blocks is at least d. Indeed, consider 
the new right end block. The outermost parts, which correspond to the lower 
d components of the 2d dimensional block vector are d dimensional since they 
correspond to the orthogonal image of the previous d-dimensional right-most 
block. The same holds for the new left-most block. Since the new extremal 
blocks are again orthogonal, their dimension is equal to d. 

This concludes an induction argument. □ 



4. Examples 

From the results presented in the previous section, it turns out that the ac- 
tual challenge in the construction of matrix filters is the choice of the orthogonal 
matrices M k in the d(2d — l)-dimensional Lie group 0(2d). 

To give examples of such construction, we restrict ourselves to the case of 
the special orthogonal group SO(2d) of rotations. 

4-1. Construction of a two-channel filter bank 

In this first example, we let d = 2 and consider a representation of 5*0(4) in 
terms of products of Givens rotation matrices G(£,m,6), with 1 < I < m < 4, 
whose only non-zero elements are given by: 

gkk = 1, k ^ £,m gu = g mm = cos 9, g lm = -g mt = sin 9. 

Let us take, as initial set of QMF filters, the Haar matrix filter bank, where: 

A(0) = A(l) = B(0) = -B(l) = J. 

We first consider a 6-parameter rotation transformation on the the following 
sequence of 2 x 2 matrices: 










1 I 



I -I 










This transformation produces two 4-length sequences A(k), B(k), k = 0, . . . , 3 
depending on 6 parameters. They can then be arranged as 








A(Q) B(0) 
A(l) B(l) 



A{2) B{2) 
A(3) B(3) 








and transformed through another 6-parameter rotation. The resulting filter 
bank consists of two 6-length sequences A(k),B(k), k = 0, . . . , 5 depending 
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on 12 parameters. These degrees of freedom can be exploited to impose addi- 
tional constraints. In particular, we can require the full rank constraint |2]) and 
the second order sum rule condition (which produces an additional vanishing 
moment on the wavelet filter): 

J2(-l) k kA(k) = 0. 

fcez 

As an example, let S(6), with 6 — [0%, #a> #3> 04, #5> be the 6-parameter 
rotation matrix obtained as the following product of Givens matrices: 

S(6) = G(l,2,9 4 )G(3,4,0 3 )G(2,3,9 2 )G(l,A,di)G(l,3,e 6 )G(2,4,9 5 ), 



explicitly given by 

C 6 CiC 4 - S 6 S 2 S4 
-C 6 ClS 4 - S 6 S 2 C4 

-s 6 c 2 c 3 - c e sis 3 

SqC 2 S 3 - C 6 SiC 3 



— S5S1C4 + C5C2S4 
S5S1S4 + C5C2C4 

-C5S2C3 - S5C1S3 
C5S2S3 - S5C1C3 



S 6 CiC 4 + C 6 S 2 S4 
— S e CiS4 + CeS 2 C4 

cqc 2 c 3 - S 6 SiS 3 
-c 6 c 2 s 3 - S 6 SiC 3 



C5S1C4 + S5C2S4 
— C5S1S4 + S5C2C4 
-S5S2C3 + C5C1S3 

s 5 s 2 s 3 + C 5 CxC 3 



with 

c-i = cos(^), Si = sin.(0»), i = 1, — ,6. 

Let now Mq = S((f>), Mi = S(ip) be the rotations (depending on the 12 free 
parameters <f>i, . . . , 4>e, ipi, . . . , tpo) respectively used for the first and the second 
transformation. 

The 12 degrees of freedom can be fully exploited to impose the 12 full 
rank/sum rule conditions. The solution is given in terms of the following pa- 
rameter vectors: 

4> = [-1.530817, -2.054355, -2.642328, 0.495166, 1.413293, 1.728299], 
V> = [-2.345058, 2.382453, -1.422064, -1.696487, 1.165227, -1.439620], 

which give rise to filters convergent to the matrix scaling function and wavelet 
illustrated in Fig. [2j 

4-2. Construction of a ^-parameter family 

In this subsection we propose a special construction of a family of filters with 
d = 2 and support length < 4. 

Instead of working with the group 50(4) itself, we consider a representa- 
tion of its elements in terms of infinitesimal generators, given by the following 
elements of its Lie algebra: 

X a = E*> j - E j '\ l<i<j<4, a = a(i, 3) 

where E*' J is the matrix whose elements are all zero except the element at the 
i-th row and j-th column, which is equal to one. 
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Figure 2: The four components of the 2-channcl scaling function (left) and wavelet (right) 
associated to the parameters obtained in the example in Section 4.1. 



We take the following ordering of the indices for the corresponding 

subscript a = a(i,j), 1 < a < 6: 

(1,2), (1,3), (1,4), (2,3), (2,4), (3,4). 

By taking exponentials of linear combinations, we obtain the elements of the 
group as 



where we have set £ = £2, • • • , S,6] T G K 6 . 

Let us now take two rotations parameterized by £ and Starting from the 
Haar system, whose polyphase matrix is given by: 



H(z) 



I I 

I -I 



our construction gives us a 12 parameter family of orthogonal filters with pa- 
i 12 . The symbols of these filters are given by 



rameters [£ T ,£' T ] T € 



zl 

1 



M(£') 







l I 



M{i)H{z). 



Not all of them are different however. And not all of them arc full rank ei- 
ther. Indeed, consider two matrices in U, V € 5*0(2). The following choice of 
matrices in SO (A) will all yield the Haar system 



M = 



For small perturbations of the Haar system a linear analysis of the space of 
filters obtained that way may be performed. For this we write 

L(t,f;z) = z^L-^t 1 ) + z°L (t,?) + (£,£')• 



u 





, M' = 


■ v T 








V 
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This defines an embedding of the space of such filters into M 3 ' 16 by taking the 
coefficients of the matrices Lk, k= — 1, 0, 1. 

Consider now the Jacobian of the nonlinear mapping (£,£') i-> £(£,£'; •) — > 
R 3 ' 16 at the origin and call this linearization T. As elementary linear algebra 
shows, the rank is 10 and hence the defect is 2. A basis of its kernel is given by: 

kernT = span{7 1 ,7 2 } 

with 

7l = [-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, If 72 = [0, 0, 0, 0, 0,-1, 1, 0, 0, 0, 0, 0] T 

So at least locally we obtain a parameterization in terms of ]R 12 /kcrnr. 

We take as complementary system the following 10 vectors written as matrix 
of column vectors 



1 
































1 
































1 
































1 
































1 





























1 
































1 
































1 
































1 
































1 
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1 






























Note that the two directions 7l and 7 2 correspond exactly to the above noted 
family of matrices which do not alter the Haar wavelets. 

We now want to explore those directions that yield one parameter subgroups 
of full rank filters. Only a sufficient condition will be given. The condition of 
full rank requires that 

£(€,€';1) = H(1). 

We are unable to solve this system in general. However, restricting our research 
to special solutions of the form L(t£,t%';z), t e R, we can construct a finite 
dimensional submanifold of full rank filters spanned by the one-dimensional 
abelian subgroup of rotations in SO (4) x 5*0(4) with their generators forming 
a sub- vector space (of its Lie algebra) . 

Since for t = the system is full rank, it is enough to require that the 
derivative with respect to t is for all t: 

±L{tt,t?;l) = 0. 

Substituting this and deriving with respect to t yields the following matrix 
system of equations: 
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I 

1 



I 

1 



\ a 



I I 




' 


" 


I -I 











I 

1 



A sufficient condition is given by the inner parenthesis to vanish. This yields 
a linear system of 16 equations for the 12 parameters £, Due to linear 
dependencies among the equations the system has a 6 dimensional solution 
space which is spanned by the following 6 vectors 
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The family of full rank filter banks that are generated by rotations associated 
with the corresponding infinitesimal rotations however is only 4 dimensional. 
Indeed, as may be verified through elementary linear algebra, the intersection 
of the linear span of these vectors and the 10 dimensional space of directions of 
local filter systems is a four dimensional vector space spanned by the following 
column vectors, which, for the convenience of the reader, we have splitted into 
the respective part for £' and £ corresponding to the outer and inner rotation 
matrices: 







i 



















-1 

















1 
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1 











1 























If we denote with , the matrices of dimension 6x4 corresponding to the 
above column vectors and let p — [t] 7 8 7 uj,(} T be the vector containing the free 
parameters, then the final form of our parameterization is expressed in terms of 
the following vectors: 



f = Gep=[ o c v w 



£ = G e p= [ -o C to n e o 
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An expression for the filters A(z) and B{z) depending on the above four pa- 
rameters is however difficult to give, due to the presence of the exponentials of 
the matrix sums J2 a €otX a , J2 a C a X a - 

Strong simplifications occur when we set all but one direction parameters to 
zero. All of the corresponding one-parameter filter families can thus be explicitly 
given. 

The family corresponding to the only non-zero direction parameter rj is given 
in terms of the following symbols: 



A(z) = (z + l) 



a 2 z + (l -a 2 ) (1 
z(l ~ z)aVl - a 2 (1 



z)ayT—a? 
a 2 )z 2 + a 2 z 



(10) 



B(z) = (l-z) 



a 2 z-(l-a 2 ) (1 + z)ay/l - 



—z(z + l)av 1 — a 2 (1 — a 2 )z 2 — a 2 z 
where we have set a = cos r) (see Fig. 3) . 

Observe that B(z) is obtained as PA(—z)P, where P - 



(11) 



1 

-1 



As to the convergence of such filters to matrix scaling function/ wavelets, ob- 
serve that the autocorrelation symbol C(z) = \A H (z)A(z) is given by C(z) = 
(z+i) j rp^ e p 0S j^j ve definiteness of C(z) assures the convergence of A(z) to 
an L 2 (M.) d orthogonal matrix scaling function as proved in 0]. 

Note that essentially the same family is attained taking u as unique non zero 
parameter. In such case only the role of £ and £' exchange. 






Figure 3: The four components of the 2-channel scaling function (left) and wavelet (right) 
associated to the symbols (jlOj, (TTTJ with a = cos(n/6) 

Let us now consider the case of a second family, obtained by letting 6 be the 
only non zero parameter in £, By setting a = cos(y26), A(z) = | (z+l)A(z), 
B(z) = 4 (1 — z)B(z) the elements of A(z) and B(z) read as: 

AiiO) = (a- l) 2 z 2 + {3 + 2a-a 2 )z 

A 12 (z) = z{l - z){l - a) (l + a - y/2(l - a 2 )) 
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A 21 (z) = (1 - z)(l - a) (l + a - Zy /2(1 - a 2 )) (12) 
^22(2) = (a + 1) (2 - 2a - ^2(1 - a 2 )) z 2 + (3a 2 + 2a - l)z 
+ (3 - a 2 - 2a + (a + 1)^2(1 - a 2 )) 



B u (z) = (a - l) 2 z 2 + (a 2 - 2a - 3)z 

B 12 (z) - z(z + 1)(1 - a) (l + a + y/2{l - a 2 )) 

Bai(«) = (z + - a) (-1 - a + z^C 1 - « 2 )) 
B 22 (z) = (a + 1) ^2 - 2a + ^2(1 - a 2 )) z 2 + (1 - 3a 2 
+ (3 - a 2 - 2a - (a + l)y/2(l - a 2 )) 



(13) 



2a)z 



Elementary computations show that also in this case the corresponding au- 
tocorrelation symbol C{z) is positive definite, which indicates the convergence 
of the filters to orthogonal matrix scaling functions/ wavelets for any value of 
the parameter a in [—1,1] (see Fig. 4). 






Figure 4: The four components of the 2-channel scaling function (left) and wavelet (right) 
associated to the symbols \12\ , ( |13| with a = 

The last one parameter family is that obtained by setting all parameters 
except C to zero. However, this is a trivial case, which means that, setting 
a = sin(C), the symbols A(z) and B(z) reduce to the diagonal forms: 



A{z) = (z + 1) 



B{z) = (l-z) 



(a 2 - ay/ 1 - a 2 )z + 1 - 2a 2 + (a 2 + ay/ 1 - a 2 );?" 1 
1 



(14) 



(a 2 + ay/ 1 - a?)z + 2a 2 - 1 + (a 2 - ay/l - a 2 );?" 1 
1 



(15) 
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and the scheme thus reduces to a combination of two scalar schemes, namely 
the Haar scalar system plus a parameterized version of the Daubechies system 
(attained for a — v3/2). This last system is convergent for any a G [—1, 1]\{0}. 

We conclude the section presenting an example of a 2 parameter family 
obtained by considering the situation # = £ = 0. In this case, the argument of 
the exponential function is a sum of two commuting matrices and thus admits 
a simple expression as product of the exponentials of each matrix. As a result, 
the explicit expressions of A(z) and B(z) = PA(—z)P in terms of the two 
parameters a = sm(r]), b = sin(a;) are: 



A(z) = (*+!) 



b 2 z 2 + (1 - a 
(1 - z) (aVT^- 



- b 2 )z + a 2 



(1 - z) (bVl - b 2 z + aVl-a 2 ) 



+ (1 - a 2 - b 2 )z + b 2 



(16) 



B(z) 



+ (1 



(1 + z) (b^/l^W - aVl-a 2 +) 



(1 + z) (aVT- 
-a 2 z 2 + (1 



by/T 
-b 2 )z 



b 2 z) 
-b 2 



(17) 



It is interesting to note that the case a = b = 1/2 produces a symbol A(z) 
with an extra (z + 1) factor, which contains the quartic B-spline symbol (which 
gives rise to a non orthogonal scalar scheme) on the diagonal. However this 
scheme is essentially diagonal, in the sense explained in |3! , and it is equivalent 
to two independent Daubechies scalar schemes. 





Figure 5: The four components of the 2-channel scaling function (left) and wavelet (right) 
associated to the two parameter symbols JTfjjl, (17) with a = 0.4, b = 0.6. 



5. Conclusions 

This paper discusses a way to parameterize families of d x d matrix wavelet 
filters of full rank type based on the one-to-one correspondence between QMF 
systems and orthogonal operators which commute with the shifts by two. It 
provides some specific examples, in particular the special construction obtained 
in terms of elements of the Lie algebra of SO(2d). Explicit expressions for the 
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filters in the case d = 2 are given, as a result of a local analysis of the parame- 
terization obtained from perturbing the Haar system. This strategy can be used 
also to generate filters with d > 2. Indeed, the idea of the construction is very 
general, and can even apply to the bivariate nonseparable setting. Furthermore, 
in all cases, the parameterization can be exploited to realize matrix systems tai- 
lored to the specific applications. Future research work will be carried on in this 
direction. 

Moreover, we are pursuing a wide experimentation connected to the appli- 
cation of multichannel wavelet filters, including those constructed in this paper, 
which shows the advantages of this tools versus traditional scalar wavelet tech- 
niques. The first results will appear in a forthcoming paper. 

Appendix A. Filter coefficients 

In this appendix, the explicit expressions of the filters taps derived in the 
Example section are given (see Tables A.1-A.4). 
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Table A.l: Coefficients of the matrix filters associated to the functions plotted in fig 



k 



A(k) 



0.054311209333498010694 
0.047827189604400026824 
0.087708930185872110064 
-0.088952962252079719437 
0.65527236381392715834 
-0.35145174509107688357 

1.088476921748056695 
-0.07789144192861210554 
0.29041642689141252135 
0.30362455490066007842 
-0.17618585197276649584 
0.16684440342858505769 



-0.16684440354507635408 
-0.17618585179038651232 
-0.30362455474026860115 
0.29041642693165752616 
0.07789144188237754498 
1.0884769216695398767 
0.35145174524818843912 
0.65527236401729840172 
0.088952962017450863579 
0.087708930001853805185 
-0.047827189524548346775 
0.054311209170036901976 



B(k) 



0.0543112093345566594 
0.0478271896021416066 
0.0877089301796203924 
-0.088952962258738763 
-0.276396211477648534 
-0.602889006833760122 

" 0.248606647088857046 
1.11448373738822482 

-0.290416426894247737 

-0.303624554889251552 
0.176185851960222128 

-0.166844403430672834 



-0.166844403526512642 
-0.176185851776288344 
-0.303624554700031024 
0.290416426927014915 

1.1144837374540959 
-0.248606647031728134 
-0.602889006560141438 " 
0.276396211518095014 
-0.0889529620232064194 
-0.087708930001465138 
0.0478271895207618912 
-0.0543112091753285287 
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Table A. 2: Coefficients of the matrix filters associated to the symbols \10) , with —1 < 



k 


A(k) 


B(k) 





1 - a 2 ay/1 — a 2 




1 - a 2 — asj\ — a 2 















1 

a\f\ — a 2 a 2 
a 2 — as/l - a 2 



i\/l — a 2 —a 2 
a 2 a\/l — a 2 



-aVl — a 2 1 — a 2 



— a 2 a 2 — 1 



Table A. 3: Coefficients of the matrix filters associated to the symbols \12) , \13\ , with —1 < 
a < 1. The choice a = leads to the functions plotted in Fig. |4| 

k AA{k) 



I -a 2 -2 a + (a + 1)^2(1 - a 2 ) +3 -a 2 
3 + 2a -a 2 (a - 1) (-1 -0+^/2(1 - a 2 ) 

(a - 1) y/2(l - a 2 ) 2 a 2 + 2 + (a + 1)^2(1- a 2 ) 



-1 + a 2 a 2 - (o+ 1)^/2(1 - a 2 ) + 1 + 2a 
(a-1) 2 -(o-l) (-l-a+ ^2(1- a 2 )) 

-(a -1)^2(1- a 2 ) -(a + 1) (-2 + 2 a + ^2(1 - a 2 )) 



-1 + a 2 -2 a- (a + l)v/ 2 (! - a 2 ) +3 -a 2 
-3-2a + a 2 - (a - 1) (l + a + ^2(1 - a 2 )) 

-(a -1)^2(1- a 2 ) -2 a 2 -2 + (a + 1)^/2(1- a 2 ) 
4 
1-a 2 a 2 + (a + 1)^2(1 - a 2 ) + l + 2a 
-(a-1) 2 (a-1) (l + a+v/2(l- a 2 )) 

(a-1) ^2(1 - a 2 ) - (a + 1) {2 - 2 a + ^2(1 - a 2 )) 
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Table A. 4: Coefficients of the matrix filters associated to the symbols \1G) , ( | 1 7| , with — 1 < 
a,b< 1. The choice = 6 = 1/2 leads to the functions plotted in Fig.p] 



A(k) 



a 2 a^l — a 2 



bVT^¥ b 2 

1 - b 2 bVT^S 2 

a\/l — a 2 1 — a 2 

1 — a 2 — a\/l — a 



-bVT^¥ 1 - & 2 
6 2 -^VT^fe 2 



— a 2 a 2 



B(k) 



— ay/1 — a 2 



Wl - fe 2 6 2 

-l + b 2 bVT^¥ 

a\/l - a 2 -1 + a 2 

1 — a 2 a\/l — a 2 
6\/l — fe 2 1 - 6 2 



-6 2 -bs/l^W 

2 



flVl — a- 
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